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Intisari 

Penelitian untuk membuat desain Electrical Capacitance Tomography dilakukan menggunakan metode nu- 
merik Finite Difference Method (FDM) untuk forward problem dan ditentukan tingkat akurasinya terhadap 
metode analitik. Solusi berbasis FDM yang diperoleh, digunakan dalam proses rekonstruksi menggunakan 
metode Landweber dan Newton-Raphson. Perbandingan citra hasil rekonstruksi dari kedua metode dianalisis 
berdasarkan resolusi dan akurasi citra yang dihasilkan. Hasil penelitian ini menunjukkan bahwa akurasi FDM 
dengan melakukan pemotongan garis pada objek domain, didapatkan 90,57% pada garis horizontal dan 87,43% 
pada garis vertikal. Adapun hasil untuk rekonstruksi didapatkan akurasi proses rekonstruksi 98% untuk Newton- 
Raphson dan 88% untuk Landweber. Resolusi citra yang direkonstruksi dengan menggunakan New ton-Raphson 
mampu membedakan 2 benda dengan jarak pisah minimal 3 cm, sedangkan Landweber mampu membedakan 
2 benda dengan jarak pisah minimal 8 cm. Hasil akurasi dan resolusi tersebut berlaku untuk sistem ECT 8 
elektroda dengan objek domain lingkaran. 


Abstract 

The Electrical Capacitance Tomography system has been designed using the Finite Difference Method as a for¬ 
warding problem-based solution. Accuracy of the solution has also determined with the analytical method. The 
reconstruction proses using Landweber and Newton-Raphson method has been conducted so that a comparison 
of resolution and accuracy of the resulting image have obtained. The results of this research, indicate that the 
accuracy of the Finite Difference Method by cutting a line on a domain object was obtained at 90.57% on the 
horizontal line and at 87.43% on the vertical line. The result of the reconstruction process based on the Newton- 
Raphson method produced an accuracy of 98%, while the Landweber method only gained 88%. Two objects 
can be distinguished well using the Newton-Raphson method when they are at least 3 cm apart. Meanwhile, the 
other method is only able to distinguish two objects, at a minimum distance of 8 cm. 
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I. PENDAHULUAN 

Tomografi adalah teknik pengamatan yang mampu 
menggambarkan sebuah objek secara cepat dan tidak bersi- 
fat merusak atau mengubah struktur benda [1]. Berbagai 
macam aplikasi dari tomografi yang telah diterapkan adalah 
menggambarkan anatomi manusia [2], mendeteksi gangguan 
tulang [3], menentukan defect pada material [4], menentukan 
kualitas makanan [5], dan monitoring regulasi dari pipa 
industri [6]. Salah satu jenis tomografi adalah Electrical 
Capacitance Tomography (ECT), yang memanfaatkan sifat 
listrik yaitu kapasitansi. ECT bekerja dengan memanfaatkan 
distribusi kapasitansi yang diukur pada elektroda di per- 
mukaan objek, untuk menentukan distribusi permitivitas 
objek yang diamati. Kelebihan tomografi jenis ECT adalah 


implementasi yang sederhana karena bekerja pada frekuensi 
rendah [7]. Agar implementasi dari ECT lebih luas, maka 
dilakukan sebuah penelitian untuk membuat sistem ECT 
dan menganalisis tingkat resolusi dan akurasi citra yang 
dihasilkan pada proses rekonstruksinya. 

Pembuatan sistem ECT dilakukan menggunakan dua 
tahapan. Tahapan pertama yaitu forward problem dengan 
menentukan distribusi potensial berdasarkan distribusi per¬ 
mitivitas, dan yang kedua adalah invers problem dengan 
menentukan distribusi permitivitas berdasarkan distribusi 
potensial [8]. Penyelesaian forward problem dilakukan 
dengan menggunakan solusi numerik agar didapatkan dis¬ 
tribusi potensial pada kondisi permitivitas tertentu. Solusi 
numerik pada penelitian ini digunakan Finite Difference 
Method (FDM) atau metode beda hingga untuk menentukan 
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Gambar 1: a). Model objek tanpa benda b). Model objek dengan 1 benda ukuran 20 x 20 cm c). Model objek dengan 2 benda 

masing-masing ukuran 10 x 10 cm dengan jarak pisah. 


distribusi potensial pada sistem ECT yang digunakan. FDM 
merupakan metode numerik untuk menentukan nilai potensial 
di suatu titik yang berdasarkan pada titik-titik disampingnya 

[9] . Kelebihan metode ini adalah dapat menentukan solusi 
masing-masing titik dengan mudah dan dapat menggam- 
barkan distribusi potensial keseluruhan yang bersifat distkrit 

[10] . Pada penelitian ini, solusi numerik yang didesain akan 
dibandingkan dengan solusi analitik, sehingga didapatkan 
tingkat akurasi dari solusi numerik. 

Bagian selanjutnya, untuk permasalahan invers problem 
atau rekonstruksi, digunakan dua metode yaitu Landweber 
dan Newton-Raphson , untuk mendapatkan hasil perbandingan 
kedua metode tersebut pada tingkat akurasi dan resolusinya. 
Metode Newton-Raphson memiliki prinsip meminimasi nilai 
kuadrat error dari selisih pengukuran kapasitansi dengan pen- 
dekatan nilai kapasitansi, sehingga didapatkan hasil rekon¬ 
struksi yang lebih akurat [11]. Sedangkan metode Landwe¬ 
ber , memiliki prinsip optimalisasi dengan selisih hasil rekon¬ 
struksi terukur dan hasil rekonstruksi melalui pedekatan, un¬ 
tuk mendapatkan hasil rekonstruksi yang akurat [12]. Dua 
metode tersebut sudah pernah dilakukan untuk mendesain 
sistem tomografi, seperti rekonstruksi ECT menggunakan 
Newton-Raphson [13], dan rekonstruksi ECT menggunakan 
Landweber [14]. Pada penelitian ini, akan dilakukan rekon¬ 
struksi menggunakan dua metode tersebut, dan selanjutnya 
akan dianalisis tentang tingkat akurasi dan resolusi yang di- 
hasilkan. 


II. METODOLOGI 

Penelitian ini dibagi menjadi beberapa tahapan yaitu 
pembuatan model objek, penyelesaian forward problem , 
penyelesaian invers problem , dan analisis akurasi dan resolusi 
citra hasil rekonstruksi. Masing-masing tahapan memiliki 
peranan penting untuk membuat sistem ECT dan proses 
analisisnya. 


Pembuatan model objek 

Pembuatan model objek ECT dilakukan dengan mendesain 
sistem ECT berbentuk lingkaran dengan 8 elektroda yang 


disusun pada permukaannya dengan sumber tegangan 50 volt, 
seperti yang ditunjukan pada Gambar 1. 

Model objek yang ditunjukan pada Gambar 1 menggu¬ 
nakan 8 elektroda, dan selanjutnya akan dibuat cacahan de¬ 
ngan jumlah grid persegi 50 x 50 cm 2 untuk mendesain so¬ 
lusi numerik dan hasil rekonstruksi pada titik-titik grid- nya. 
Model objek dengan permitivitas yang dijadikan background 
adalah 1, sedangkan benda memiliki permitivitas 7, dengan 
ukuran 20 x 20 cm 2 dan 10x10 cm 2 pada jarak pisah 10 cm. 
Pengukuran dilakukan untuk satu elektroda sebagai sumber 
sedangkan 7 elektroda lain digunakan untuk menentukan kap¬ 
asitansi. Proses tersebut dilakukan dengan pergantian sumber 
untuk elektroda yang lain. 

Penyelesaian forward problem 

Penyelesaian forward problem digunakan untuk menentukan 
distribusi potensial f yang diperoleh dengan menyelesaikan 
persamaan Poisson [15]. 


V • cV ( h = —p 


( 1 ) 


dengan p adalah rapat muatan dan e adalah permitivitas benda. 
Selanjutnya, persamaan 1 dikorelasikan dengan variabel ka¬ 
pasitansi sehingga mendapatkan persamaan kapasitansi yang 
bergantung pada distribusi permitivitas e, untuk dua elektroda 
yang memiliki beda potensial A Vij [16]. 

Cij = --A— J^e(x,y)V^(x,y)(U (2) 

Nilai kapasitansi Cij merupakan nilai yang diukur pada se- 
tiap pasangan elektroda, dimana Cij merupakan fungsi in¬ 
tegral permitivitas dan gradien potensial untuk koordinat x 
dan y terhadap d £, sehingga penelitian ini menghasilkan 56 
hasil pengukuran kapasitansi untuk model objek yang digu¬ 
nakan. Berdasarkan Pers.(2), dapat diambil sebuah matrik 
yang menghubungkan antara variabel kapasitansi dan variabel 
permitivitas, yang disebut sebagai matrik sensitivitas Sij [17]. 


-L 


V4>i • V4>. 
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a) b) 


Gambar 2: a). Distribusi potensial dengan solusi analitik b). Distribusi potensial dengan solusi numerik FDM. 


Penyelesaian invers problem atau rekonstruksi 

Terdapat dua metode yang digunakan pada proses rekon- 
straksi yaitu metode Newton-Raphson dan metode Landwe- 
ber. Persamaan yang digunakan dalam metode Newton- 
Raphson ditunjukkan oleh Pers.(4) [18]. 

Gk+1 =G k + (ST •Sk + jiy 1 Si (<Ck - C) (4) 

Sedangkan persamaan dari metode Landweber ditunjukkan 
oleh Pers.(5) [19]. 

G k+ i = G k + a k »S T (C-S*G K ) (5) 

Kedua metode tersebut digunakan untuk menentukan dis¬ 
tribusi permitivitas hasil rekonstruksi G berdasarkan hasil 
pengukuran nilai kapasitansi C. Kedua metode akan diber- 
lakukan pada jumlah iterasi 1 hingga 100.000 untuk mendap- 
atkan hasil rekonstruksi yang baik. 

Analisis akurasi forward problem , akurasi proses 
rekonstruksi, dan resolusi 

Akurasi solusi numerik yang diperoleh menggunakan FDM 
diperbandingkan dengan solusi analitik. Seleksi terhadap 
model objek yang akan ditentukan nilai akurasinya, dilakukan 
dengan cara mengambil daerah hasil slicing di bagian garis 
tengah objek domain. Hasil slicing akan menghasilkan garis 
horisontal dan garis vertikal untuk ditentukan nilai akurasinya. 
Perhitungan error dari solusi numerik mengunakan FDM di¬ 
lakukan dengan Pers.(6) dengan V adalah variabel potensial. 

error = |FmtmeHfc ~ x i 00 % (6) 

^analitik 

Adapun error untuk proses rekonstruksi ditunjukkan pada 
Pers.(7). 

error = x 100% (7) 

lie'll 2 

Berdasarkan data error di atas, tingkat akurasi ditentukan 
dengan menggunakan Pers.(8). 


Pengukuran resolusi ditentukan dengan mengamati jarak 
pisah minimum 2 buah benda hasil rekonstruksi yang masih 
dapat dibedakan satu sama lain. Jarak pisah diambil mulai 
dari jarak 10 cm dan diperkecil dengan interval 1 cm. Jarak 
minimum yang diperoleh menunjukkan kemampuan sistem 
ECT untuk membedakan 2 benda yang terpisah. 

III. HASIL DAN PEMBAHASAN 

Hasil penelitian ini dibagi menjadi 3 bagian. Bagian per- 
tama yaitu hasil penyelesaian forward problem , bagian kedua 
yaitu penyelesaian invers problem , dan bagian ketiga adalah 
resolusi citra hasil rekonstruksi. 

Hasil penyelesaian forward problem 

Penyelesaian forward problem menggunakan FDM dan solusi 
analitik pada kasus sederhana atau permitivitas homogen 
berupa distribusi potensial yang ditunjukkan pada Gambar 
2. Solusi analitik pada Gambar 2 didapatkan dari penyele¬ 
saian persamaan Laplace dengan memberikan syarat batas 
v(x,0) = 50 sin(7rx/50) volt, v(x,oc) = 0, dan ekspansi deret 
menggunakan deret Fourier [19]. Sedangkan solusi numerik 
menggunakan syarat batas v(x,0) = 50 sin(7rx/50) volt. Syarat 
batas tersebut dipilih karena fungsi sinus akan menghasilkan 
potensial tertinggi pada posisi tengah, dimana posisi tersebut 
adalah posisi dari elektroda pertama pada sistem ECT yang 
dibuat. 

Hasil akurasi solusi numerik terhadap solusi analitik pada 
titik-titik yang ditentukan, memiliki akurasi rata-rata yaitu 
90,57% pada garis horisontal, dan 87,43% pada garis vertikal. 
Hasil akurasi FDM yang diperoleh dapat dikatakan tinggi 
karena memiliki tingkat akurasi lebih dari 80%, sehingga so¬ 
lusi numerik FDM dapat digunakan untuk mendesain matrik 
sensitivitas pada sistem ECT. 

Hasil penyelesaian invers problem 

Penyelesaian invers problem menggunakan dua metode yaitu 
Newton-raphson dan Landweber. Hasil rekonstruksi dua 
metode tersebut ditunjukkan pada Gambar 3. 


akurasi = 100 % — error 


( 8 ) 
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Gambar 4: a) Benda asli berjarak pisah 8 cm. b) Hasil rekonstruksi Newton-Raphson berjarak pisah 8 cm. c) Hasil rekonstruksi 
Landweber berjarak pisah 8 cm. d) Benda asli berjarak pisah 7 cm. e) Hasil rekonstruksi Newton-Raphson berjarak pisah 7 cm 
. f) Hasil rekonstruksi Landweber berjarak pisah 7 cm. g) Benda asli berjarak pisah 3 cm. h) Hasil rekonstruksi 
Newton-Raphson berjarak jarak 3 cm. i) Hasil rekonstruksi Landweber berjarak pisah 3 cm. j) Benda asli jarak 2 cm. k) Hasil 
rekonstruksi Newton-Raphson jarak 2 cm. 1) Hasil rekonstruksi Landweber jarak 2 cm. 
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Hasil pada Gambar 3 menunjukkan bahwa metode Newton- 
Raphson dapat merekonstruksi citra lebih baik daripada 
Landweber , dengan tingkat akurasi pada Newton-Raphson 
mencapai 98%, sedangkan untuk Landweber mencapai 88%. 
Selain tingkat akurasi, berdasarkan Gambar 3 dapat dilihat 
bahwa rekonstruksi Newton-Raphson memiliki bentuk yang 
lebih mirip benda asli daripada rekonstruksi Landweber. 
Namun, kedua metode memiliki hasil rekonstruksi dengan 
efek smoothing. Efek ini menghasilkan sebuah citra yang 
memiliki perbedaan kontur yang memiliki arti perubahan 
nilai permitivitas pada permukaan benda, sehingga efek 
tersebut akan mengganggu hasil rekonstruksi yang dihasilkan 
oleh sistem ECT. 

Selain efek smoothing yang dihasilkan, citra hasil rekon¬ 
struksi kedua metode memiliki kelemahan berupa keadaan 
tidak dapat menggambarkan benda lancip atau memiliki 
sudut. Hal ini dapat dilihat pada citra hasil rekonstruksi 
bagian ujung benda, dimana kedua metode tidak dapat 
menggambarkan bentuk benda secara sempurna. Selain 
itu, hasil rekonstruksi Landweber masih menghasilkan citra 
yang berbentuk mirip lingkaran, hal ini berbeda seperti yang 
dihasilkan oleh Newton-Raphson , dimana Newton-Raphson 
hanya membentuk lengkungan pada ujung sisi dari benda 
yang diamati. Adapun hasil dari efek smoothing dan kelema¬ 
han dalam merekontruksi benda yang memiliki sudut, citra 
yang dihasilkan dengan permitivitas 7 mengalami reduksi 
ukuran. Benda asli dengan ukuran 20 x 20 cm 1 2 3 4 5 6 7 pada 
permitivitas 7 hanya dapat direkonstruksi menjadi benda 


berukuran 12 x 12 cm 2 untuk permitivitas 7. Sedangkan 8 
x 8 cm 2 sisanya merupakan hasil dari efek smoothing yang 
diperoleh pada proses rekonstruksi. 

Analisis resolusi 

Resolusi yang didapatkan merupakan tingkat kemampuan 
proses rekonstruksi membedakan dua benda yang terpisah de¬ 
ngan jarak tertentu. Hasil citra dua metode yang digunakan 
untuk menentukan resolusi dapat dilihat pada Gambar 4. 

Hasil pada Gambar 4 menunjukkan bahwa metode Newton- 
Raphson dapat membedakan dua benda dengan jarak > 3 cm, 
sedangkan pada metode Landweber dapat membedakan dua 
benda dengan jarak > 8 cm. Oleh karena itu, resolusi yang 
dihasilkan oleh metode Newton-Raphson lebih baik daripada 
yang dihasilkan oleh metode Landweber pada penelitian ini. 

IV. SIMPULAN 

Berdasarkan hasil penelitian yang telah dilakukan, dida¬ 
patkan bahwa solusi numerik menggunakan FDM dapat digu¬ 
nakan untuk menyusun sistem ECT karena memiliki akurasi 
> 80% terhadap solusi analitik pada kasus permitivitas 
homogen. Hasil rekonstruksi citra menggunakan metode 
Newton-Raphson memiliki akurasi dan resolusi yang lebih 
tinggi dibandingkan dengan metode Landweber , untuk sistem 
ECT yang didesain pada penelitian ini. 
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